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SUMMARY 


A  numerical  treatment  of  the  influence  of  the  aeroelastic  deformations 
on  the  longitudinal  stability  and  control  characteristics  of  an  aircraft 
with  a  swept-back  wing  of  large  aspect  ratio  is  carried  out  for  the  sub- 
critical  speed  range.  Only  the  deformation  of  the  wing  and  the  bending 
of  the  fuselage  have  been  taken  into  consideration. 

The  method  used  to  calculate  the  deformations  of  the  wing  is  based  on 
the  superposition  of  a  number  of  pre-selected  load  distributions  (A 
Superposition  Method  for  Calculating  the  Aeroelastic  Behavior  of  Swept 
Wings,  Aug. 1951.  by  R.B.  Brown  et  alii).  An  investigation  of  the 
accuracy  of  the  method  has  been  made  by  comparing  the  numerical  results 
with  those  of  another  more  exact  method  which  proved  the  reliability  of 
the  present  method. 


As  a  result  of  the  study  it  was  found  that  the  stability  and  control 
of  the  aircraft  are  only  slightly  affected  by  the  aeroelastic  deformation. 


SOMMAIRE 

Un  traitement  num&ique  de  1'  influence  des  deformations  adro^lastiques 
sur  la  stability  longitudinale  et  sur  les  caract6ristiques  de  contrdle 
d’un  avion  &  voilure  en  fleche,  d’ allongement  4lev^,  est  effectu4  en  ce 
qui  concerne  la  gamme  des  vitesses  sous  critiques.  Seules  la  deformation 
de  la  voilure  et  la  flexion  du  fuselage  ont  4t4  prises  en  consideration. 
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Le  mode  de  calcul  des  deformations  alaires  utilise  se  base  sur  la 
superposition  d’un  certain  nombre  de  repartitions  des  charges  predetermines 
(A  Superposition  Method  for  Calculating  the  Aeroelastic  Behavior  of  Swept 
Wings,  Aug. 1951,  by  R.B.  Brown  et  alii).  La  precision  de  ce  mode  a  et 4 
verifies  en  confrontant  les  resultats  nuaeriques  obtenus  et  ceux  obtenus 
d’aprds  un  autre  node  plus  precis,  ce  qui  a  perais  de  mettre  en  evidence 
la  sflrete  du  mode  actuelleaent  employe. 


Par  suite  de  cette  etude  il  a  M  constate  que  la  stabillte  et  le 
contrftle  de  1*  avion  ne  sont  que  peu  sensibles  4  la  deformation  aero- 
eiastique. 
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THE  INFLUENCE  OF  AEROELASTICITY  ON  THE  LONGITUDINAL 
STABILITY  OF  A  SWEPT  WING  SUBSONIC 
TRANSPORT  AIRCRAFT 

C.M.  Kalkman* 


1.  INTRODUCTION 

By  order  of  the  “Netherlands  Aircraft  Development  Board' ,  a  quantitative  calcula* 
tion  was  carried  out  for  a  hypothetical  aircraft  with  elastic  swept-back  wings  in 
order  to  investigate  the  changes  in  the  longitudinal  stability  and  control  character¬ 
istics  as  a  consequence  of  the  deformations.  Only  the  influences  of  wing  deflection 
and  fuselage  bending  are  taken  into  account. 

The  determination  of  the  load  distribution  on  the  elastic  wing  under  consideration 
has  been  performed  by  the  method  of  superposition.  Many  other  methods  of  calculation 
may  be  found  in  the  literature  on  this  subject;  most  are  complicated  and  time- 
consuming.  The  method  of  superposition  is  relatively  simple,  so  that  it  seemed 
desirable  to  investigate  its  reliability  from  the  point  of  view  of  accuracy.  The 
influence  of  compressibility  of  the  air  was  introduced  into  the  calculation  by  an 
approximate  expression.  To  avoid  too  much  numerical  calculation  the  reliability  of 
the  method  of  superposition  was  checked  at  only  one  Mach  number  (M  =  0.8). 

With  the  help  of  the  method  of  superposition  the  distributions  of  lift  and 
pitching  moment  along  the  elastic  wing  are  determined  for  various  values  of  M  and 
the  dynamic  pressure  (horizontal  flight  at  two  altitudes).  With  these  data  the 
quantities  which  are  of  importance  for  the  static  longitudinal  stability  and  control, 
e.g.  location  of  the  centre  of  pressure  and  lift  gradient,  are  determined.  More¬ 
over,  the  downwash  at  the  horizontal  tail  is  calculated  by  using  the  distribution 
of  lift  along  the  span  of  the  elastic  wing. 

With  the  parameters  obtained  in  this  way,  elevator  angles  to  trim  in  level  flight, 
elevator  deflections  per  ‘g’ ,  and  the  corresponding  manoeuvre  margins  were  calculated 
as  functions  of  speed.  To  be  able  to  detect  the  aeroelastic  effects  these  calcula¬ 
tions  were  carried  out  for  the  rigid  aircraft  and  the  flexible  aircraft  with  and 
without  taking  into  account  the  influence  of  compressibility  of  the  air. 


2.  STATEMENT  OF  THE  PROBLEM 

The  purpose  of  this  investigation  is  to  determine  to  what  extent  the  static 
longitudinal  stability  and  control  characteristics  of  an  aircraft  with  slender  swept- 
back  wings  are  affected  by  aeroelastic  deformations.  The  behaviour  of  large  trans¬ 
ports  and  bombers  was  of  primary  interest  in  this  respect  as  large  effects  are  to  be 
expected  in  these  cases. 


*  Nationaal  Luchtvaart laboratorium,  Amsterdam ,  Netherlands 
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To  get  an  impression  of  the  quality  and  magnitude  of  the  above-mentioned  factors 
quantitative  calculations  were  carried  out  for  a  hypothetical  aircraft,  the  shape  and 
dimensions  of  which  roughly  agree  with  those  of  the  Boeing  B-47.  In  order  to  limit 
the  scope  of  the  calculations  to  a  reasonable  extent  only  wing  deflections  and  aero- 
elastic  fuselage  bending  are  considered. 

•  For  this  purpose  the  aircraft  was  supposed  to  consist  of  the  following  parts: 

(a)  An  elastic  wing  without  twist  with  a  symmetrical  cross-section  and  an  angle 
of  sweep  of  the  quarter-chord  line  6  -  35°  ,  an  aspect  ratio  k  =  8  and  a 
taper  ratio  of  1/3.  The  corresponding  wing  area  amounts  to  162  m2; 

(b)  A  flexible  fuselage  of  about  20  m  length; 

(c)  A  rigid  all-movable  tail  with  -  35°  and  k  =  4  (area  =  25  m2). 

The  separate  effects  of  deformation  may  be  compared  by  combining  a  rigid  wing  with 
a  flexible  fuselage  and  a  flexible  wing  with  a  rigid  fuselage.  Moreover,  the 
influence  of  compressibility  has  been  taken  into  account;  but  only  sub-critical 
Mach  numbers  have  been  considered. 

A  general  view  of  the  problems,  which  present  themselves  successively  when 
calculating  the  aeroelastic  influence  on  longitudinal  stability,  may  be  obtained 
most  effectively  with  the  help  of  the  pitching-moment  equation  for  the  equilibrium 
condition.  This  equation  can  be  approximately  expressed  as  follows: 


where 


v 


t 


v 


a.  c. 


1 


t  At  St 
A  S 


(2) 


and 


v 


a.  c. 


v 


c.g. 


(3) 


The  lift  gradient  A  of  the  wing  and  of  the  horizontal  tail  At  ,  the  location 
of  the  wing  centre  of  pressure  xa> c< /c  and  the  downwash  gradient  de/da  are 
considered  to  be  affected  by  aeroeiastic  deformation  in  Equation  (1).  If  it  is 
possible  to  calculate  these  quantities  for  different  values  of  Macb  number  and 
dynamic  pressure  the  effects  of  deformation  of  wing  and  fuselage  can  be  determined 
immediately. 

A  discussion  of  the  applied  calculation  methods  will  be  given  in  the  following 
sections. 
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3.  METHODS  TO  DETERMINE  THE  AERODYNAMIC  LOADS  ON  ELASTIC  WINGS 

3. 1  General 

The  difficulty  which  presents  itself  when  calculating  aerodynamic  loads  on  lifting 
surfaces  is  a  consequence  of  the  mutual  influence  of  load  and  deformation.  Most 
existing  methods  have  the  disadvantage  that  complicated  iteration  procedures  are 
necessary  to  determine  the  equilibrium  condition  of  the  structure.  Moreover,  these 
methods  are  often  based  on  ‘strip  theory’ ,  so  that  induced  aerodynamic  effects  are 
not  taken  into  account. 

In  Reference  1  a  method  is  given  -  the  so-called  ‘superposition  method’  -  which 
leads  rapidly  to  the  required  results.  Since  the  accuracy  of  this  method  is  some¬ 
what  adversely  affected  by  the  way,  in  which  the  influence  of  compressibility  is 
taken  into  account,  a  second  more  exact  method  is  applied.  By  comparing  the 
numerical  results  of  these  two  methods  the  feasibility  of  the  superposition  method 
can  be  judged. 

3.2  The  Superposition  Method 

3.2.1  Description  of  the  Method 

Although  it  is  rather  difficult  to  calculate  the  final  angle-of-attack  distribution 
for  a  wing  with  a  given  initial  (unloaded)  a-distribution,  it  is  quite  possible  to 
calculate  the  angle-of-attack  distribution  for  the  unloaded  condition  if  the  final 
angle-of-attack  distribution  (or  final  loading)  is  given.  This  reverse  procedure  is 
the  basis  of  the  superposition  method. 

Starting  from  a  number  of  arbitrarily  selected  final  angle-of-attack  distributions 
(ctp)  ,  the  load  distributions  for  the  wing,  considered  to  be  rigid  in  this  condition, 
are  now  calculated.  The  deformations  (tXg)  which  are  a  consequence  of  these 
calculated  loadings  may  be  determined-,  if  the  stiffness  parameters  of  the  wing  are 

known. 


The  angle-of-attack  distribution  dj  of  the  unloaded  wing  then  follows  from: 

dp  =  dj  +  dE  (4) 

This  relation  applies  to  any  wing  section.  Prom  each  selected  distribution  of  the 
final  angle  of  attack  (dp)  follows  a  corresponding  one  along  the  unloaded  wing. 

A  number  of  dp-distributions  are  considered  now  from  which  the  corresponding 
dj-distributions  are  calculated.  The  given  dj-distribution  may  then  be  approximated 
by  the  following  expression: 

dj  =  A^ij  +  Bgdj^  +  Cgaj^  + .  (5) 

In  general  the  following  four  dp-distributions  are  chosen: 
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(a)  A  distribution  equal  to  the  given  aj-distribution; 

(b)  Three  additional  distributions  ap^,  a?  and  Op3>  proportional  respectively 
to  -rj  ,  7j2  and  rj3  (17  =  (2y/b)cos£  ,  the  dimensionless  coordinate). 


The  resulting  aerodynamic  loadings  cause  the  elastic  deformations  of  the  angle- 
of-attack  distributions  a&  and  agi  (i  =  l  to  3)  ,  which  may  be  approximated  by 
the  following  cubic  equations  in  17  : 


«E  =  aPfaji?  +  b^2  +  CjT)3) 


Og  =  qP(aTj  +  b7j  +  C7j3) 


V  (6) 


These  equations  show  moreover  that  Og  is  linearly  dependent  on  the  dynamic 
pressure  q  and  the  quantity  P  ,  which  represents  the  ratio  between  the  lift 
gradients  with  and  without  compressibility  effects.  The  relation  between  ctg  and 
P  is  based  on  the  assumption  that  the  load  distribution  does  not  change  shape  with 
increasing  speed  below  the  critical  speed.  This  assumption,  which  is  made  in 
Reference  1  without  motivation,  appeared  to  be  acceptable  according  to  comparative 
numerical  results. 


Furthermore,  if  we  suppose  that  Hooke’ s  law  applies  and  that  the  aerodynamic 
loading  is  linearly  dependent  on  the  angle  of  attack,  the  following  relations  may 
be  derived  from  Equations  (4)  and  (5): 


“E  =  +  +  cs“e3  + 

=  Vrt  +  Vf,  +  cb°t3  + 


(7) 


As  As  ,  Bg  ,  CB .  follow  from  the  approximation  given  fay  Equation  (5)  the 

final  angle-of-attack  distribution  can  now  be  determined  by  starting  from  a  given 
dj-dlstribution;  the  resulting  aerodynamic  properties  may  then  be  derived. 


The  given  az -distribution  is  now  composed,  according  to  Equation  (5),  in  the 
following  way  of  the  aj-distrlbutlons  corresponding  with  (a^r  +  a^): 

AsaI  +  BsaI  +  csaI  +  aj  =  aj  (8) 

12  9 


where 


“f  =  aI 


(») 


Vlth  the  help  of  Equations  (4)  and  (9),  Equation  (8)  nay  be  transformed  into 

Vlj  +  B.al2  +csaI,  =  “« 


(10) 
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With  this  relation  the  'superposition  coefficients'  nay  be  calculated  along  the 
lines  given  in  the  next  section. 

Once  the  coefficients  Ag  ,  Bg  and  Cfl  are  known,  the  quantities  which  are 
of  importance  for  this  investigation,  viz.  wing  lift  gradient,  wing  pitching  moment 
and  downwash  gradient  at  the  tall,  may  be  calculated  from  the  following  equations: 


The  first  terms  on  the  right-hand  sides  of  these  equations  represent  the  coefficients 
of  lift  gradient,  pitching  moment  and  downwash  gradient  respectively  at  the  horizontal 
tail  for  the  rigid  wing.  The  second,  third  and  fourth  terms  refer  to  the  corresponding 
quantities  occurring  as  a  consequence  of  the  selected  a-distributlana  along  the  rigid 
wing  span  (a  =  77,  17s  and  173)  and  may  be  considered  as  correction  terms. 

3. 2.2  The  Calculation  of  the  Superposition  Coefficients  (Ag,  Bg  and  Cg) 

The  superposition  coefficients  A,  ,  Bb  and  C,  must  be  derived  from  Equation  (10) 
as  mentioned  in  the  preceding  section.  To  get  an  expression  for  the  aj-distributions 
the  corresponding  deformations  agt  should  be  calculated  from  the  loading  caused  by 
the  selected  distributions  ari  .  In  agreement  with  the  fact  that  these  a^-dlstributions 
are  of  the  third  degree  at  most  (ctr1  =  rj,  r?  and  tj3)  the  number  of  superposition 
coefficients  does  not  amount  to  more  than  three-  The  elastic  angle-of-attack 
distributions  og4  have  been  approximated  by  cubic  equations  in  17  : 

a*i  =  qP(ai7j  +  b^3  +  c^3)  (12) 


The  product  qP  is  dependant  on  flight  altitude  and  Mach  number.  For  this 
investigation  two  altitudes  (1000  n  and  10,000  n)  and  a  Mach  number  range  up  to  0.8 
have  been  considered  (sub -critical  speeds). 


The  influence  of  compressibility  is  expressed  hy  the  quantity  P  ,  the  ratio  of 
the  lift  gradients  approximated  according  to  Reference  S  for  rigid  wings  in  oon- 
pressible  and  Incompressible  flow: 


(Ar)g  _  h  +  2  corn# _ 

(Ar)fc0  ~  \/(l  -  M*  cos*  0)  +  2  coo  6 


(IS) 


The  aagle-of-attaok  distributions  aXi  can  now  be  calculated  with  the  help  ef 
Equation  (4): 


Sfj  =  +  “ft 
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Inserting  the  above-sent loned  expressions  for  and  api  as  functions  of  T) 

leads  to: 

aIt  =  V  “  <jbi  =  (1  *  »!«?)■>?  -  bjQPrj*  -  c^Pi?3 


al2  =  732  •  at2  ~  *  a2*>Pj?  +  (1  -  baqP)7j2  -  c2qPrj3  j-  (14) 


a!  =  t}3  -  “g  =  -  a3qPr?  -  b3qPTj3  +  (1  -  c3qP)i73 

3  3  J 

The  quantity  derived  fros  (=  aT)  should  be  approxisated  by  a  cubic  in 
7)  as  sell  as  the  quantities  : 

a{  =  a-rj  +  bv)2  +  C7)3  (15) 

Equating  (10) 

Vi,  +  B,ai,  +  c.ai,  =  «*t 

12  3 

and  Equation  (15)  for  any  value  of  17  leads  to: 


Pros  those  equations  the  unknown  quantities  A,  ,  Bg  and  Cg  oan  be  determined  for 
each  combination  of  altitude  and  Mach  number. 


4.  CALCULATION  OP  THE  AEBOPYNAHIC  LOAOINO  ON  TEE  BIOID  UNO 

TO  determine  the  deforaatloos  (a,  )  it  is  necessary  to  knot  the  loadings  caused 
by  the  seleoted  distributions  of  the  final  angle  of  attack  afl  along  the  span  of 
the  wing  (considered  to  be  infinitely  stiff).  There  are  different  methods  available 
to  calculate  these  lift  distributions.  The  method  developed  by  Multhopp*  has  been 
used  for  the  present  Investigation.  An  Important  advantage  of  this  method  ooegared 
with  most  others  of  the  mm  reliability  from  the  point  of  view  of  aeenraoy  is  the 
relatively  small  extent  of  the  required  oempmtatioaal  work.  This  is,  maoag  other 
things,  qua  to  a  oorreet  seleetiom  of  the  location  of  the  eo-eelled  pivotal  potato 
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on  the  lifting  surface.  In  this  respect  the  number  of  coefficients  of  the  summation 
terms  representing  the  downwaah  integral  is  reduced  considerably,  if  the  pivotal 
points  along  the  span  are  selected  in  such  a  way  as  to  form  the  projections  of  the 
points  obtained  by  dividing  the  half  circle  on  the  wing  span  into  equal  parts. 

The  number  of  pivotal  points  (m)  applied  in  this  investigation  amounted  to  25, 
according  to  Multhopp’ s  recommendation  (wing  aspect  ratio  X  -  8): 

m  >  3\/l  -  Mz  (17) 

The  method  results  in  a  number  of  linear  equations  with  the  local  lift  and  the  local 
pitching  moment  as  unknowns. 

Apart  from  the  smaller  amount  of  computational  work  this  method  also  has  the 
advantage  that  it  applies  to  the  swept-back  wing  and  that  it  is  accessible  to 
compressibility  calculations.  When  applying  the  superposition  method  use  was  made 
exclusively  of  the  results  for  the  incompressible  airflow  condition,  while  the 
'exact*  method  carried  out  for  comparative  purposes  (  at  M  =  0. 8)  was  based  on 
Multhopp* s  method  for  the  compressible  case. 

5.  CALCULATION  OF  THE  DEFOBMATIONS  DUE  TO  A  6IVEN 

AERODYNAMIC  LOADING 

Computation  of  the  aerodynamic  loading  mentioned  in  the  preceding  section 
provides  values  for  the  dimensionless  lift  per  unit  length  of  wing  span  and  for  the 
pitching  moment  at  the  pivotal  points.  With  these  data  the  resulting  deformations 
may  be  calculated.  An  oblique  co-ordinate  system  was  chosen  for  this  purpose,  being 
most  suitable  for  the  swept-back  wing  configuration.  The  X-axis  coincides  with  the 
axis  of  symmetry  of  the  wing  and  the  Y-axis  with  the  elastic  axis,  which  is  supposed 
to  be  a  straight  line. 

The  deformation  of  the  wing  may  be  given  hy  the  following  differential  equations: 

5“  =  9m"x  + 

(17) 

^—2  =  Q21M,  +  ^22*** 

The  angles  <px  and  qy  represent  the  rotations  around  the  X-axis  and  the  Y-axis 
respectively  In  an  arbitrary  section  y  of  the  wing.  They  are  components  along  the 
X-  and  Y-axis  of  the  resulting  rotation  7  •  while  Mx  and  My  are  the  normal 
projectioas  on  those  axes  of  the  resulting  moment  vector  acting  at  this  section. 
Accordingly  the  coefficients  Qtl  in  Equation  (17)  can  be  expressed  as  follows: 
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1 

Qn  =  Bcos20 


Q 


22 


1  1  2 
-  +  -  tan  A 
T  B 


(18) 


1  tan  A 

Q  j ,  =  Q21  =  •  —  - — 

1 2-  21  B  cos  8 

The  quantities  B  and  T  represent  the  bending  and  torsional  stiffnesses 
respectively  of  a  section  normal  to  the  elastic  axis,  of  which  the  angle  of  sweep 
is  equal  to  6  .  To  make  a  numerical  evaluation  of  the  deformations  possible, 
certain  assumptions  referring  to  the  stiffnesses  of  the  wing  had  to  be  made.  Prom 
a  rough  estimation  of  the  bending  stiffness  (B0  =  EI0)  at  the  wing  root  followed 
Bfl  =  10 12  kgcm2  .  The  desired  data  were  partly  derived  from  publications  on  the 
Boeing  B-47  aircraft  and  partly  obtained  from  estimations  on  a  theoretical  base. 


In  order  to  realise  a  reasonably  perceptible  deformation  of  the  lifting  surface 
it  was  supposed  that  the  bending  stiffness  along  the  span  of  the  tapered  wing 
varied  proportionally  to  the  fifth  degree  of  the  chord  (taper  ratio  of  the  wing  is 
equal  to  1/3).  Furthermore,  the  assumption  was  made  that  for  each  wing  section  the 
ratio  of  bending  stiffness  to  torsional  stiffness  (T)  is  constant  (a  practical 
value  of  B/T  =  4  has  been  chosen). 


The  first  equation  in  (17)  is  related  to  the  bending  of  the  lifting  surface 
(around  the  X-axis)  while  the  second  equation  refers  to  the  change  in  the  angle 
of  attack.  Since  only  the  latter  is  of  importance  with  respect  to  this  investiga¬ 
tion  attention  will  be  mainly  confined  to  this  equation. 


The  moment  Mx  in  the  section  y  may  be  represented  by  the  following  relation: 


«x 


/ 


(b/2)/cos  8 
ky ^ (y j  -  y)  cos  8  dyx 


(y  <  yt) 


(19) 


where  kyt  is  the  loading  per  unit  length  (Y-axis)  at  the  section  yt  .  The 
expression  under  the  Integral  sign  gives  the  moment  of  the  loading  of  a  wing  strip 
yt  around  an  axis  in  the  section  y  parallel  to  the  X-axis. 

Accordingly  the  moment  My  in  the  section  y  may  be  expressed  by: 


"y 


/ 


(b/2)/coa  9 


(20) 


where  m7l  is  the  moment  around  the  Y-axls  per  unit  length  along  this  axis  in  the 
section  .  The  moments  following  from  calculations  according  to  Multhopp  are 
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referred  to  an  axis  through  the  quarter-chord  point  of  the  section  and  normal  to 
tjie  plane  of  symmetry  of  the  wing.  Consequently,  these  moments  should  be  reduced 
to  that  point  of  the  section  where  it  iB  intersected  by  the  elastic  axis;  further¬ 
more  they  should  be  multiplied  by  cos  8  before  they  can  be  substituted  in 
Equation  (20). 

Substituting  Equations  (19)  and  (20)  in  the  second  Equation  in  (17),  while  the 
co-ordinates  y  and  y^  are  made  dimensionless  [(2ycos#)/b  =  77]  ,  the  following 
expression  may  be  developed: 

Bcp  -b2tan£  f  b  (\  1  ,  \  f 

-  -  /  k  (77  -  77)d77  +  - -  -  +  -  tan 2  8)  m_  dr).  (77  <  77  ) 

3y  4B  cos 2  8  /  vi  1  1  2cos0\T  B  //  1 

J -n  x  Jn 


Suppose  that,  according  to  the  above,  the  bending  and  torsional  stiffnesses 
vary  along  the  span  as  the  fifth  degree  of  the  chord,  and  taking  a  taper  ratio  of 
1/3  gives,  we  obtain: 


B  =  B.(l  -  2/3  77) ! 


T  =  T.(l  -  2/3  77)' 


Integration  of  Equation  (21)  yields: 


b3tan$ 


“o'*8  ^ 


k  (Vi  -  '»7>dT71 


(1  -  2/3  77)  ■ 


b  /I  1 

+ - *—  —  +  —  tan2  8 

4cosz0\T  B 


0  l  Jv 


drj 

B_  d T) ,  - - - 1 

17 »  1  (i  -  2/3  77)- 


This  equation  can  be  integrated  by  using  influence  functions.  For  the  sake  of 
brevity  only  the  final  result  is  presented  here: 

0,(77)  =  9  (77)  cos  8  J 


,  ,  77  b*qtan0  V 

to*w)>‘  =  '  55  i^sZ  a,“jyj  * 
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The  Influence  functions  i21  and  i22  are  expressed  in  the  spanwise  co-ordinate 
\p  (17  =  cos  yp)  to  allow  immediate  multiplication  with  Multhopp’s  dimensionless 
coefficients  for  local  lift  and  pitching  moment  y  and  y.  .  The  relation  between 
these  coefficients  is  given  by: 


=  2qbyj 

(■q  )j  =  2qbc(/x  +  0. ly)j  cos  8 


(25) 


For  an  arbitrary  section  parallel  to  the  plane  of  symmetry  the  chord  is  given  by: 

(26) 


3J  =  co^  -  \c0*^ 


Finally  the  influence  functions  occurring  in  Equation  (24)  can  be  written  as 
follows:  ] 

i81('/'.'/'1>  =  3/4 


1 

.  -  2 

9(1  -  2/3  cosi/O 

1 

1 

-  1  - 

(1  -  2/3  cost/Zj)3 

16 

(1  -  2/3  COS*/:.)1* 

1  « 

*••**•»'*  3/»{  (1  .  2/300.^.  -1} 


(27) 


These  equations  apply  to  the  case  for  which  *p  >  \pl  ,  if  </<  <  the  quantity  i px 
in  the  above  expressions  for  the  influence  functions  has  to  be  substituted  by  'p  . 


The  summation  in  Equation  (24)  has  to  be  taken  over  all  Multhopp  stations  on  one 
half  of  the  wing  (n  =  13)  . 


For  the  sake  of  completeness  the  vertical  wing  displacements  are  derived  from  the 
first  Equation  in  (17).  To  get  an  insight  into  the  degree  of  deformation,  wing  deflec¬ 
tions  were  calculated  only  at  ths  tip  (for  which  1  =  13)  from: 


ztip  =  zi»is 


e  »» 

w  qb8  ^ 

=  W2B oCC^X  (i‘l>UrJ 

W  J«1 


77  qb*c0tand 

96  B.cos*0  *  l,^iJ^J 


(28) 


6.  CALCULATION  OF  TEE  DOVNWASH  OF  THE  ELASTIC  VINO 
AT  THE  HORIZONTAL  TAIL 

Because  the  superposition  method  is  based  on  ths  superposition  of  a  number  of 
quantities  corresponding  with  different  angle-of-attack  distributions,  it  will  be 
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necessary  to  calculate  the  downwash  caused  by  each  selected  basic  load  condition 
(angle-of-attack  distribution)  in  order  to  be  able  to  determine  the  resulting  down- 
wash  at  the  horizontal  tail.  The  station  for  which  the  downwash  was  determined  is 
located  at  an  equal  distance  from  the  plane  of  symmetry  as  the  first  pivotal  (Multhopp) 
point  on  the  wing.  This  location  produces,  according  to  practice,  a  representative 
value  of  the  average  downwash  at  the  horizontal  tail  (roughly  1/6  of  the  tailspan 
out  of  the  centre),  so  that  the  computational  work  can  be  limited  to  this  station. 

The  downwash  at  this  point  due  to. a  given  a-distribution  may  be  considered  to 
consist  of  4  parts: 

<  =  el  +  e2  +  e3  +  (29) 

The  components  £  t  and  e2  are  due  to  the  bound  vortices  of  the  right  and  left 
wings  respectively,  while  £3  and  £„  originate  from  the  trailing  vortices  of  the 
right  and  left  wings  respectively. 

The  4  quantities  £t  ...  eH  are  calculated  with  the  help  of  the  following  formula 
for  the  induced  normal  velocity  w  in  a  point  at  a  distance  r  from  a  vortex-line 
with  circulation  P  : 


r  ds  cos  a  r  cos  a 

dw  =  - s—  =  - —  da  (30) 

4nr  47 rl 

where  a  is  the  angle  between  r  and  the  normal  I  on  the  vortex  line. 

Integration  of  Equation  (30)  yields  formulae  for  the  components  £t  . 

For  the  downwash  angle  £t  induced  hy  the  bound  vortex  of  the  half  wing  located 
at  the  sane  side  of  the  plane  of  symmetry  as  the  chosen  tail  station  P  the 
following  relation  was  derived: 

e  _  b2*co»2g  *  V  _ Vjito'Pj 

1  877  24  /([Gib cost/',  -  (Isin0  +  p)}*  +  l*coa*0]3) 

i«i  1 
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The  symbols  l  and  p  are  defined  in  the  sketch  herewith.  The  distance  of 
the  station  P  from  the  bound  vortex  at  the  quarter-chord  line  of  the  wing  is 
denoted  by  l  ,  while  the  distance  of  this  point  P  to  the  plane  of  symmetry  is 
equal  to  p  . 

The  downwash  angle  e2  induced  by  the  bound  vortex  on  the  half  wing  at  the 
other  side  of  the  plane  of  symmetry  may  be  expressed  by: 


b^l'cos2#  77  ^7 


£2  = 


_ 7j  B*n0i 

/([{ftbcoscp,  -  (Vain  8  -  p)}z  +  l^cos2#])^ 


877  24 


(32) 


ial 


l<Pi 


Evidently  the  quantity  V  represents  the  distance  of  P  to  the  quarter-chord  line 
of  the  other  half  wing. 

In  the  sketch  an  arbitrary  trailing  vortex  is  drawn. 

The  relation  for  the  downwash  angle  e3  Induced  at  the  point  P  by  the  trailing 
vortices  at  the  right  hand  side  of  the  plane  of  symmetry  reads: 


«.-b2 


Ay_\ 

Wh 


1  + 


L  -  -cost p.  tan 6 
2  1 


‘  Icos^i(Lt“5  +  p)  +  +  L*) 


3y 


b 

— sint/t^ 


—cost p.  -  p 
2  1 


•i 


(33) 


The  differential  quotient  By/ 3 y'  can  be  derived  directly  from  the  Multhopp 
computations  (  3y/3y'  =  Jib  dy/drj)  . 

The  downwash  angle  f ,  Induced  by  the  trailing  vortices  at  the  other  side  of  the 
wing  can  be  written: 


L  •  —  cost/tj  tan  8 
2 


*  =  »«^3y ^  +  bc0B^i(p : htm0i  7 p, ;  L.) 


b 

-sln^ 


M 


a  , 
—cost Pi  4-  P 
2 


(34) 


The  suanatioo  £  is. carried  out  over  the  Multhopp  pivotal  points  selected  before. 


For  each  of  the  selected  angle-of-attaok  distributions  one  value  of 
e  =  («t  +  cf  +  «,  +  «,)  is  oomputed.  It  is  possible  now  to  detersine  a  final  value 
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of  the  downwash  gradient  with  the  superposition  method  starting  from  a  chosen  initial 
condition. 

The  following  aircraft  dimensions  have  been  used  for  the  numerical  evaluation  of 
the  downwash  formulae: 

L  =  19  m 
b  =  36  m 
8  =  32.5° 

/  b  tt\ 

p  =  2.35  m  =  -sin  — 

V  2  24 J 


7.  FUSELAGE  BENDING 

In  addition  to  the  aeroelastlc  deformation  of  the  wing,  fuselage-bending  has  also 
been  taken  into  account.  For  this  purpose  use  was  made  of  the  approximative  method 
given  in  Reference  5,  in  which  a  relation  is  derived  between  the  lift-gradients  of 
the  horizontal  tailplane  with  and  without  the  effects  of  elastic  deformation  of  the 
fuselage: 


A 


t 


AtR 

1  +  AtRJq 


(35) 


Here  At  and  AtR  represent  the  lift  gradients  of  the  horizontal  tail  with  and 
without  the  influence  of  fuselage  bending.  The  quantity  J  can  be  roughly  estimated 
with  the  help  of  the  following  relation: 

S  l 2 

J  =  0.715  — * -  (36) 

<EIf>0 

in  which  it  has  been  assumed  that  the  bending  stiffness  of  the  fuselage  varies 
according  to: 


The  X-axis  in  this  case  was  supposed  to  coincide  with  the  axis  of  symmetry  of  the 
fuselage.  The  length  of  the  fuselage  lt  (=  16.5m)  is  taken  from  the  origin  of  the 
X-axis,  where  the  fuselage  is  assumed  to  be  clamped  (roughly  at  the  quarter-chord 
point  of  the  wing  nidchord).  The  height  of  the  fuselage  varies  parabolical ly. 
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The  numerical  value  of  the  Quantity  J  is  8.65  x  10* *  kg" Si2,  based  on  the  follow* 
ing  estimation  of  the  bending  stiffness: 

(Elf) 0  =  6  x  10*  kgm2 

With  these  data  the  lift  gradient  At  can  now  be  calculated  for  each  Mach  number. 

8.  THE  ‘EXACT*  METHOD  TO  CALCULATE  THE  DISTRIBUTION  OP  THE 

AERODYNAMIC  LOADING  OF  AN  ELASTIC  WING 

In  the  preceding  sections  it  was  mentioned  that  in  order  to  check  the  results  of 
the  superposition  method  a  second  method  has  been  applied  to  calculate  the  load 
distribution  on  the  elastic  wing  more  exactly. 

This  method  is  based  on  the  solution  of  a  number  of  linear  equations  with  the 
dimensionless  lift*  and  pitching-moment  coefficients  y  and  /x  as  unknowns.  The 
coefficients  of  y  and  /x  in  these  equations  are  obtained  hy  using  the  following 
fundamental  equation: 


Op  (Vi)  =  (38) 

In  this  equation  the  angle*of -attack  distribution  Ogfy)  is  a  linear  polynomial 
of  y  and  /x  .  The  final  load  distribution  on  the  elastic  wing  is  determined  by 
the  final  angle-of*attack  distribution  Oy  ,  by  means  of  the  Multhopp  equations 
for  fx  and  y  .  which  can  be  written  for  a  wing  with  a  symmetrical  croas*sectioo 
as  follows: 

ni.  nii 

7l  =  an(l{  -  +^>8^  +2Clj/i3  (3#> 

Jml  j»l 

In  the  same  way  an  expression  is  obtained  for  /x  : 

"il  % 

fh  =  *ii<Bi  *  ■i'>°*i  +2  EljMj  <40> 

j»i  j»i 

The  subscripts  on  the  summation  symbol  are  introduced  to  denote  that  i  t  j  . 

All  coefficients  occurring  in  Equations  (39)  and  (40)  can  be  calculated  from  special 
tables  given  by  Multhopp,  when  the  geometry  of  the  wing  is  known.  Subsequently 
the  final  angle*of*attack  distribution  ay  is  substituted  according  to  Equation  (38) 
by  the  Initial  and  the  elastic  angle-of*attack  distributions.  For  the  Og-distributioo 
the  expression  from  Equation  (24)  can  be  substituted,  while  the  ax*  distribution  is  a 
known  function  of  the  spsnwlse  co-ordinate.  (In  this  investigation  ax  is  constant 
along  the  span).  After  carrying  out  the  above  mentioned  substitutions  there  appear 
a  number  of  equations  which  are  linear  in  /j.  and  y  .  Solution  of  these  equations 
yields  the  final  distributions  of  lift  and  pitching  moment. 
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Although  the  principles  of  this  method  are  relatively  simple,  the  introduction 
of  compressibility  when  considering  different  air  speeds  produces  a  considerable 
quantity  of  additional  computational  work. 


9.  COMPARISON  OP  THE  RESULTS  OP  WING  LOAD  CALCULATIONS  KITH 

THE  SUPERPOSITION  METHOD  AND  THE  ‘EXACT*  METHOD 

The  results  of  the  lift  computations  obtained  with  the  ‘exact*  method  were  used 
as  a  measure  for  the  reliability  of  the  superposition  method,  dealt  with  in 
Section  3.2. 

The  most  important  influence  of  the  differences  between  the  results  of  both 
methods  is  due  to  the  way  in  which  compressibility  is  taken  into  account.  The 
superposition  method  is  based  on  the  assumption  that  the  shape  of  the  load 
distribution  does  not  change  with  increasing  Mach  number.  One  of  the  consequences 
of  this  assumption  is  that  the  centre  of  pressure  of  the  section  does  not  shift. 

The  exact  method  applies  transformed  equations  for  the  computation  of  lift  and 
pitching  moment  according  to  the  *Prandtl-Glauert*  rule  when  compressibility  is 
taken  into  account.  The  effect  of  this  transformation,  the  so-called  'co-ordinate- 
reduction*  ,  is  a  change  of  shape  in  the  lift  distribution,  especially  along  the 
wing  chord,  so  that  movements  of  the  centre  of  pressure  occur. 

The  comparison  between  the  two  methods  was  carried  out  for  a  Mach  number  of  M  =  0.8 
(sub-critical)  in  order  to  obtain  the  largest  possible  difference  between  the  results. 
As  a  matter  of  fact  this  Mach  number  is  relatively  high  as  far  as  the  validity  of 
the  ‘Prandtl-Glauert*  rule  is  concerned. 

In  Figure  1  load  distributions  are  given  for  an  altitude  of  10,000  a,  expressed 
by  the  product  2b  7  =  CjC  as  a  function  of  the  spanwise  co-ordinate  (17) 

Integration  of  this  product  over  the  span  yields: 


Since  the  calculations  discussed  here  start  from  an  angle  of  attack  of  1  radian 
the  lift  coefficient  obtained  with  Equation  (41)  is  in  fact  a  lift  gradient  A  . 

Two  curves  have  been  drawn  in  Figure  1  for  the  rigid  wing.  One  of  these  load 
distributions  (full  line)  was  obtained  by  multiplying  the  corresponding  results 
computed  for  incompressible  flow  with  the  factor  P  earlier  mentioned  (Bqn. 13), 
whereas  the  second  me  (dotted  line)  was  determined  directly  for  the  required  Mach 
number  according  to  Multhopp's  method. 

The  load  distributions  calculated  for  the  aeroelastically  deformed  wing  calculated 
with  the  superposition  method  (full  line)  and  the  ‘exact*  method  (dotted  line)  are 
given  by  the  two  lower  curves.  The  difference  between  the  results  of  both  methods 
appears  to  be  very  small  and  nay  be  practically  neglected,  notwithstanding  the  fact 
that  the  superposition  method  is  based  00  a  very  rough  estimation  of  the  compres¬ 
sibility  effects  compared  with  the  exact  method  (see  Section  3.2). 
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Further  information  referring  to  the  reliability  of  the  superposition  method  may 
be  obtained  by  comparing  the  lift  gradients  computed  from  the  load  distributions 
with  Equation  (41).  According  to  the  exact  method  the  lift  gradient  A  changes 
from  5.51  for  the  rigid  wing  to  4.73  in  the  deformed  condition.  The  corresponding 
figures  for  the  superposition  method  are  5.43  and  4.71.  On  account  of  these  results 
it  appears  that  the  applicability  of  the  superposition  method  is  permissible. 

In  a  similar  way  the  corresponding  curves  for  the  distributions  of  the  pitching 
moment  (expressed  as  2b/x  )  were  computed  and  these  are  drawn  in  Figure  2.  As  a 
consequence  of  the  way  in  which  compressibility  has  been  taken  into  account  the  two 
curves  of  the  pitching-moment  distribution  for  the  rigid  wing  show  remarkable 
differences.  However,  to  appreciate  the  influence  of  aeroelasticity  on  the  longitu¬ 
dinal  stability  of  aircraft,  it  is  the  total  shift  of  the  centre  of  pressure,  being 
identical  to  the  aerodynamic  centre  for  the  wing  with  symmetrical  cross-section  under 
consideration,  which  is  of  primary  importance,  not  its  absolute  location.  Comparing 
now  the  shifts  of  the  centre  of  pressure  according  to  both  methods,  the  following 
results  may  be  enlightening. 

According  to  the  ‘exact’  method  the  location  of  the  centre  of  pressure  for  the 
rigid  condition  is  at  -1.53  c  (that  is  behind  the  apex  of  the  wing);  this  shifts 
to  -1.49  c  for  the  elastic  wing  under  the  same  conditions.  Calculations  according 
to  the  superposition  method  produce  a  shift  from  -1.49  c  to  -1.45  c  .  In  both 
cases  a  forward  shift  is  found  of  4%  c  .  From  this  point  of  view  also  the  applica¬ 
bility  of  the  superposition  method  is  justified. 

To  judge  the  mutual  differences  between  both  methods,  the  elastic  angle-of-attack 
distribution  (a^)  has  been  considered.  The  results  are  given  in  Figure  3.  The 
total  elastic  angle  of  attack  was  computed  according  to  the  superposition  method  using 
the  following  relation: 


«E,  =  Vs,  +BsaB1  +CsaE1  +aE,  (42> 

In  this  equation  the  quantities  agi(1>  2  3)  and  represent  the  changes  of 

the  angle-of-attack  distributions  due  to  aeroelastic  deformation  at  the  Hulthopp’ 
stations  on  the  wing  caused  by  the  loadings  following  from  the  final  distributions 
“p  =  ’l.7?*.7?3  *nd  1  radian. 

The  component  curves  (dotted  lines)  of  the  resultant  ag^  (small  circles)  are 
also  drawn  in  Figure  3  in  agreement  with  Equation  (42).  To  get  an  Insight  into  the 
accuracy  of  the  approximation  of  the  elastic  angle-of-attack  distributions  Og  by 
a  third-degree  function  in  tj  (determination  of  the  coefficients  a^  ,  bj  ,  and  c,  ), 
the  resultant  etg  was  computed  moreover  from: 

Og  =  At7)  +  Bg7)2  +  Cg7)3  (43) 

The  a^-curve  is  given  by  the  full  line  in  Figure  3.  It  appears  froa  the  location 
of  the  small  circles  (Bqu. (42))  with  regard  to  this  full-line  curve  (Bqu. (43))  that 
the  approximation  is  fully  justified  as  both  computations  lead  to  the  same  result. 
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Finally,  the  elastic  angle-of-attack  distribution  according  to  the  exact  method 
was  calculated  by  substituting  the  corresponding  values  of  y  and  /z  into 
Equation  (39).  The  distribution  is  now  solved  and  the  distribution  follows 
from: 


aE1  =  aFi  *  \  <  44> 

The  values  agi  computed  in  this  way  are  inserted  into  Figure  3  as  separate  data. 
Comparing  now  the  results  of  the  two  methods  under  consideration  it  appears  from 
Figure  3  that  the  differences  are  only  small.  In  the  most  unfavourable  case  they 
only  amount  to  about  10%  of  the  total  change.  It  may  be  concluded  from  the  above 
results  that  no  predominant  objections  against  the  application  of  the  superposition 
method  are  present.  Consequently,  the  remaining  part  of  the  investigation  has  been 
carried  out  according  to  this  method. 


10.  INVESTIGATION  OF  THE  INFLUENCE  OF  AEROELASTICITY  ON  THE 
STATIC  LONGITUDINAL  STABILITY  OF  AN  AIRCRAFT  «ITH 
SLENDER  SWEPT-BACK  1INGS 

10. 1  Basic  Data  for  the  Calculation 


Before  passing  on  to  the  proper  treatment  of  the  problem  a  summary  is  given  of  the 
numerical  values  of  the  quantities  which  have  served  as  a  basis  for  the  calculations: 


W  =  100.000  kg 

5  =  162  m2 
b  =  36  m 

\  =  8  (aspect  ratio) 

taper  ratio  1/3 

cQ  =  6.75  m  (root  chord) 

6  -  31.5°  (angle  of  sweep  of 

elastic  axis) 

Table  I  gives  the  flight  conditions  for 


B0  =  10*  kgm2 
T0  =  0.25  x  10*  kgm2 

St  =  25  m2 
bt  =  10  m 

lt  =  16.5  a  (distance  of  c.g.  to  a.c. 

of  horizontal  tailplane) 

(Atr)  =  3.65 

r  M»0 

Vt  =  0.565 

which  the  deformations  have  been  calculated. 


The  centre  of  gravity  was  assumed  to  be  located  at  0.1  c  behind  the  aerodynamic 
centre  of  the  rigid  wing,  that  is  1.595  5  behind  the  apex  of  the  swept  wing.  The 
quantity  (xc>  -  x^ c< )  ,  which  occurs  in  the  pitching-moment  equation  for  the 
equilibrium  condition,'  is  for  the  compressible  as  well  as  for  the  incompressible  flow 
equal  to  0.1  5  (the  superposition  method  supposes  that  the  load  distribution  along 
the  wing  chord  does  not  change,  shape  with  increasing  Mach  number).  With  the  help  of 
the  superposition  method  the  final  distributions  of  y  and  n  were  calculated. 

Proa  these  the  pitching  moment  coefficient  (with  reference  to  the  normal  on  the  plane 
of  symmetry  in  the  apex  of  the  wing)  and  the  lift  coefficient  were  derived.  The 
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location  of  the  wing  centre  of  pressure  (xQ  =  *atC>)  was  derived  from  the  ratio 
C„/CL  .  The  numerical  values  (xc>g>  *  xcp')/c  are  given  in  Table  II. 

The  downwash  angle  at  the  horizontal  tail  also  has  been  determined  by  the  super¬ 
position  method.  The  values  de/d a  are  summarized  in  Tables  Ilia  and  Illb. 

Finally,  a  fuselage  pitching-moment  coefficient  has  been  assumed,  viz.  Cgf  =  -0.015; 
the  Influence  of  elastic  deformation  was  left  out  of  consideration  but  the  influence 
of  compressibility  was  taken  into  account  in  the  following  way: 


0.015 

✓<i  -  mV 


(45) 


10.2  Survey  of  the  Calculations 

As  has  been  mentioned  already  when  defining  the  problem, a  hypothetical  aircraft 
has  been  considered  consisting  of  an  elastic  swept-back  wing,  an  elastic  fuselage 
and  a  rigid  all-movable  tailplane. 

Because  stability  margins  have  no  practical  meaning  at  high  speed,  only  elevator 
angles  to  trim  have  been  calculated  for  the  steady  rectilinear  flight,  whereas 
manoeuvre  margins  and  elevator  deflections  per  *g'  have  been  computed  for  the  steady 
accelerated  flight. 

In  order  to, get  an  insight  into  the  separate  influences  of  aeroelasticity  and 
compressibility,  the  calculations  were  carried  out  for  the  cases  listed  in  Table  IV, 
for  two  flight  altitudes. 


10.3  Results 

The  elevator  deflections  in  steady  flight  have  been  determined  from  the  pitching 
moment  equation  for  the  equilibrium  condition.  Remembering  that  xc  p>  =  x|  C(  , 
this  equation  is,  for  the  aircraft  with  an  all-movable  tailplane, 

CM  =  0 


=  C„  +  c,  -  v*  J 


c 

In  this  equation  the  tail  volume  vt  follows  from: 

v 


M  -- 

A  \  da 


") 


)CL  +  At8t 


yt  = 


c.p. 


A  S  \  da/ 


1  + 


(46) 


(47) 


in  which 
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V  —  n 

c.p.  c.  g. 


x„  ,  -  x„  _ 

+  _ c.p. 


(48) 


The  san oeuvre  margin  may  be  written  as  follows: 

H  _  *c.b.  '  xcLg,  +  A  v  I i  (l  -  — ^  — 

*  *  c  *  fcl  A  \  da/  2^ 

and  the  elevator  deflection  per  ‘g’  may  be  calculated  from: 


(49) 


S 


g 


A*v 


ft 


(50) 


A  picture  of  the  static  longitudinal  stability  in  unaccelerated  flight  my  be 
given  with  the  help  of  the  elevator  deflections  to  trim  as  they  follow  from 
Equation  (46).  The  change  of  these  deflections  with  the  lift-coefficient  of  the 
aircraft  is  reproduced  in  Figure  4  for  Cgf  =  0  and  in  Figure  5  for 
Cgf  =  0.015// (1  -  M2)  . 

It  follows  from  these  figures  that  the  influence  of  compressibility  is  consider¬ 
ably  larger  than  that  of  aeroelasticity.  It  appears,  furthermore,  that  the  trim  change 
in  consequence  of  compressibility  alone  is  almost  constant  for  the  whole  CL-range 
considered,  while  it  is  larger  at  higher  altitudes.  This  is  in  agreement  with  the 
approximation  that  the  influence  of  compressibility  is  proportional  to  the  product 
of  CL  and  M2  (Ref.  6). 

To  get  a  better  insight  into  the  magnitude  of  the  trim  change  in  consequence  of 
aeroelasticity,  the  parameters  occurring  in  Equation  (46)  which  are  of  Interest  in 
this  respect,  viz.  (xc>gt  -  xc<p  )/5  *(=  3cL/9a)  ,  are  drawn  separately  as 

functions  of  CL  in  Figures  6  to  8. 

It  may  be  deduced  from  these  Figures  and  Equation  (46)  that  the  changes  of  the 
c.p.  location  and  of  lift  gradient  due  to  aeroelasticity  counteract  each  other. 
Considering  for  instance  in  Figure  7  the  change  of  lift  gradient  A  when  passing 
from  the  rigid-compressible  case  to  the  elastic-compressible  one,  a  large  decrease 
of  A  appears  to  occur,  which  is  largest  at  the  highest  value  of  M  .  At  the  same 
time  a  shift  in  forward  direction  of  the  centre  of  pressure  takes  place.  In 
Figure  6  a  shift  of  about  0.12  5  is  seen  at  M  =  0.8  .  The  consequence  of  these 
counteracting  effects  is  the  small  trim  change  occurring  when  changing  from  elastic 
to  rigid  conditions  in  compressible  flow.  The  same  applies  however  for  the  incompres¬ 
sible  case,  as  may  be  also  concluded  from  the  corresponding  curves. 


Physically  this  may  be  explained  also  in  the  following  way.  As  a  consequence  of 
deformation  the  c.p.  of  the  wing  will  shift  forward  causing  a  nose-up  change  of 
pitching  moment.’  Furthermore  a  considerable  decrease  of  lift  gradient  takes  place. 
TO  maintain  a  constant  value  of  CL  an  increase  of  angle-of-attack  at  the  wing  is 
required.  It  follows  now  that,  because  the  tail  angle-of-attack  increases  too,  a 
nose-down  pitching  moment  is  generated. 
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The  above-mentioned  change  of  lift  gradient  and  c.p.  shift  are  numerically  the 
most  important  parameters  with  respect  to  the  trim  changes.  Of  course  the  change  of 
downwash  angle  due  to  aeroelasticity  too  plays  a  part  in  this  respect,  but  this  is 
only  of  small  importance  compared  with  the  effects  mentioned  earlier. 

Let  us  consider  fuselage  bending  results  in  a  change  of  tailplane  lift  gradient. 

This  change  is  also  relatively  small,  so  that  only  a  small  positive  trim  change 
nay  be  noticed  from  the  curves  of  Figures  4  and  5,  when  passing  from  the  case  of  a 
flexible  wing  only  to  that  of  both  flexible  wing  and  flexible  fuselage. 

Introduction  of  a  zero-lift  pitching  moment  of  the  fuselage  (Cgf  =  -0.015^(1  *  M2) 
results  in  a  parallel  displacement  of  the  trim  curves  (compare  Fig. 4  and  Fig.  5)  for 
incompressible  flow.  An  almost  negligible  trim  change  in  a  destabilizing  sense  occurs 
for  the  compressible  flow  condition  at  high  speed.  This  effect  is  only  noticeable 
for  the  curves  referring  to  a  flight  altitude  of  1000  m. 

The  manoeuvre  margin  stick-fixed  (H^)  has  been  calculated  for  the  different  cases 
considered  using  Equation  (49)  and  drawn  as  a  function  of  CL  (or  M  )  in  Figure  9. 
Transition  from  the  rigid-incompressible  case  to  the  rigid-compressible  one  results  in 
a  decrease  of  Ha  with  II  .  The  effect  of  aeroelastic  deformation  is  larger  at 
lower  altitudes  at  the  same  Mach  number  in  consequence  of  the  larger  value  of  the 
dynamic  pressure  q  .  The  influence  of  wing  deformation  appears  to  be  only  small. 

The  manoeuvre  margin  decreases  a  small  amount  with  increasing  M  for  both  compressible 
and  incompressible  flow.  This  decrease  appears  to  be  much  larger,  however,  when  the 
fuselage  bending  is  considered  additionally  (decreasing  At  ). 

The  elevator  deflections  per  *g*  (8  )  as  a  function  of  CL  are  computed  from 

Equation  (50)  and  are  drawn  in  Figure  10  for  all  cases  considered.  The  effect  of 
compressibility  alone  is  to  cause  a  constant  shift  of  8  over  the  whole  C^-range 
with  reference  to  the  incompressible  condition  just  as  it  has  been  shown  already  for 
the  elevator  trim  curves.  This  shift  increases  with  increasing  altitude,  which  is  in 
agreement  with  an  approximate  relation  derived  in  Reference  6,  according  to  which 
this  shift  should  be  proportional  to  C^M2  . 

The  influence  of  aeroelasticity  on  Sg  ,  however,  appears  to  be  much  smaller. 

The  combined  effect  of  both  phenomena  is  relatively  large  at  low  values  of  CL  . 

At  CL  =  0.15  (altitude  1000  a)  the  elevator  deflection  per  V  is  almost  halved. 

This  effect  decreases  relatively  with  decreasing  speed. 

To  get  an  impression  of  the  magnitude  of  the  deformation  occurring  during  the 
steady-flight  conditions  considered  the  corresponding  vertical  displacements  of  the 
wing  tip  were  calculated  from  Equation  (28)  and  plotted  in  Pigure  11  together  with 
the  corresponding  ^-distributions.  It  appears  from  this  figure  that  the  displace¬ 
ments  of  the  wing-tip  (Ztlp  as  a  function  of  CL)  are  of  the  same  order  of  magni¬ 
tude  for  all  values  of  CL  ,  at  both  altitudes.  Differences  occur  only  in  the 
change  of  these  displacements  with  CL  .  The  fact,  however,  that  the  wing  lift  for 
all  cases  considered  here  must  be  nearly  equal  to  the  aircraft  weight  is  the  reason 
why  the  range  of  tip  displacements  is  only  small.  The  bending  of  the  wing  at 
10,000  a  altitude  decreases  and  at  1000  m  increases  with  increasing  CL  .  This 
is  caused  fay  the  shape  of  the  corresponding  lift  distributions  on  the  wing.  At  an 
altitude  of  1000  a  the  lift  of  the  wing  (see  Pig. 11)  is  concentrating  towards  the 
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wing  centre  with  increasing  Mach  number,  so  that  the  corresponding  bending  moment  and 
consequently  the  tip  displacement  will  decrease.  The  reverse  effect  occurs  at 
10,000  m.  It  appears  from  these  calculations  that  considerable  aeroelastic  deforma* 
tions  may  take  place  (vertical  tip  displacement  of  about  3.0  m)  without  deteriorating 
the  longitudinal  stability. 


11.  CONCLUSIONS 

In  order  to  investigate  the  possible  variations  of  the  longitudinal  stability  and 
control  of  an  aircraft  as  affected  by  aeroelastic  deformation  (particularly  of  the 
wing),  computations  have  been  carried  out  for  a  hypothetical  aircraft  with  slender 
swept -back  wings.  Of  course  the  results  only  apply  to  the  aircraft  type  considered. 

To  be  able  to  calculate  the  required  basic  data  for  this  investigation,  use  was  made 
of  a  superposition  method,  derived  by  Brown,  Holtby  and  Martin1.  With  this  method 
numerical  results  were  obtained  in  a  relatively  rapid  way  for  lift  gradients, 
c.p. -locations  and  downwash  gradients  of  the  wing  at  different  combinations  of 
dynamic  pressure  and  Mach  number. 

This  method  was  checked  from  the  point  of  view  of  accuracy  by  comparing  the  results 
with  those  of  another,  more  exact  method.  For  a  complete  application  the  latter 
requires  much  more  computational  work.  The  comparison  has  been  carried  out  for  one 
Mach  number  only  (M  =0.8)  . 

After  having  shown  the  feasibility  of  the  superposition  method  in  this  way  the 
remaining  basic  data  have  been  calculated  according  to  the  simplified  method.  With 
these  data  the  characteristics  of  the  static  longitudinal  stability  of  the  aircraft 
have  been  derived,  viz.  elevator  angles  to  trim,  elevator  angles  per  *g'  and 
manoeuvre  margins  stick-fixed.  Prom  these  final  results  the  following  conclusions 
nay  be  drawn: 

(1)  When  the  load  distributions  are  calculated  both  with  the  supeposition  method 
and  with  the  'exact  method',  it  appears  that  the  mutual  differences  of  the 
results  are  sufficiently  small  to  justify  the  application  of  the  superposition 
method  for  rapid  computations  (Fig.  1).  The  wing  lift  gradient  at  M  =  0.8  and 
at  a  dynamic  pressure  q  =  1200  kg/m2  changes  from  5.51  to  4.73  due  to  aero- 
elasticity  according  to  the  exact  method,  while  the  figures  according  to  the 
superposition  method  are  5.43  and  4.71  respectively. 

(2)  The  mutual  differences  of  the  pitching-moment  distributions,  however,  are 
much  larger  (Fig. 2).  Nevertheless  the  derived  c.p. -shifts  appear  to  be  almost 
equal  in  both  cases.  As  only  the  latter  is  of  interest  for  this  Investigation 
the  suitability  of  the  applied  method  has  been  proved  in  this  respect  as  well. 
According  to  the  exact  method  the  wing  deformation  causes  a  shift  of  the  c.p. 
from  1.530  5  to  1.490  c  ,  behind  the  apex  of  the  wing,  while  the  correspond¬ 
ing  values  according  to  the  superposition  method  are  1.495  and  1.453  5  .  In 
both  oases  a  forward  shift  was  obtained  of  about  4f  5  . 

(3)  The  deformations  of  the  wing  determined  by  both  methods  as  given  by  the  angle- 
of-attaek  distributions  also  show  small  mutual  dlfferenoss  (Fig.  3)  for 
the  flight  conditions  considered. 
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(4)  The  longitudinal  trim  changes  as  a  consequence  of  wing  deformation  of  the 
hypothetical  aircraft  considered  appeared  to  be  small  in  compressible  as  well 
as  in  incompressible  flow.  The  influence  of  compressibility  is  considerably 
larger  (Figs. 4  and  5).  The  effect  of  wing  deformation  on  longitudinal 
stability  results  in  two  opposite  effects,  viz.  a  forward  shift  of  the  c.p.  and 
a  decrease  of  lift  gradient  with  increasing  speed,  effects  which  almost 
neutralize  each  other. 

(5)  The  manoeuvre  margin  stick-fixed  only  shows  an  appreciable  decrease,  when 
fuselage  bending  is  taken  into  account  (Fig. 9).  The  influence  of  compres¬ 
sibility,  however,  is  much  larger. 

(6)  The  elevator  travel  per  ‘g*  undergoes  a  relatively  large  decrease  as  a  con¬ 
sequence  of  the  combined  influences  of  compressibility  and  aeroelasticity, 
in  particular,  at  low  altitudes  and  high  Mach  numbers  (Fig. 10).  Again  the 
influence  of  compressibility  is  dominating. 

(7)  Hie  vertical  displacements  of  the  tip  in  steady  flight  appear  to  be  of  the 
order  of  3.0  m;  only  small  variations  occur  with  changing  flight  conditions 
(see  Fig.  11). 

Notwithstanding  the  considerable  deformations  of  the  wing,  the  changes  of 
longitudinal  stability  due  to  aeroelasticity  appear  to  be  unexpectedly  small.  When 
comparing  the  above  conclusion  with  those  obtained  by  other  investigators  in  this 
field,  e.  g.  by  Campion  and  Skoog1*1 7,  similar  results  are  given.  The  counteracting 
effects  of  c.p. -shift  and  wing-lift  change  due  to  aeroelasticity  with  respect  to 
stability  and  trim  change  are  particularly  mentioned  in  both  references.  Likewise 
caution  has  been  recommended  when  generalizing  these  results  by  postulating  the  rule 
that  flexibility  of  swept-wing  aircraft  does  not  affect  stability  in  a  serious  way. 

By  a  proper  choice  of  the  stiffness  parameters  dependent  on  the  sweep  angle  of  the 
wing  it  is,  however,  possible  to  avoid  a  too  large  change  of  the  stability  character¬ 
istics. 

Contrary  to  those  publications  this  investigation  is  chiefly  confined  to  the 
aeroelastic  deformation  of  the  wing  and  fuselage  bending.  The  weight  of  the  aircraft 
wing,  which  generally  has  an  alleviating  effect  on  longitudinal  stability  changes, 
haa  been  left  out  of  consideration. 
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TABLE  I 


Plight  Conditions  Used  for  Deformation  Calculations 
(see  Sec. 10. 1) 


Altitude 

10,000  D 

1000  m 

M 

q(kg/m2) 

CL 

q(kg/m2) 

CL 

0.4 

• 

- 

1026 

0.601 

0.6 

682 

0.905 

2309 

0.267 

0.7 

927 

0.666 

3143 

0.196 

0.8 

1210 

0.510 

4105 

0.150 

TABLE  II 

Numerical  Values  of  (x.  _  -  x_  _  )/c 

§1  V •  p* 


H 

Altitude 

10,000  m 

1000  m 

mm 

Km 

BE m 

ESI 

msM 

0.8 

0.142 

wm 

TABLE  III 

Values  of  de/da  for  the  Rigid  and  the  Elastic  Ring 

(a)  Compressible  Flow 


Altitude 

M 

10,000  a 

1000  n 

Rigid 

Elastic 

Rigid 

Elastic 

0.4 

0.6 

0.7 

0.8 

0.311 

0.328 

0.352 

0.323 

0.345 

0.376 

0.290 

0.311 

0.328 

0.352 

0.305 

0.343 

0.371 

0.407 

(b)  Incompressible  Flow 


Alt  i 

tude 

10,000  m 

1000  m 

M 

Rigid 

Elastic 

Rigid 

Elastic 

0.4 

• 

• 

0.276 

0.289 

0.6 

0.276 

0.285 

0.276 

0.302 

0.7 

0.276 

0.288 

0.276 

0.308 

0.8 

0.276 

0.291 

0.276 

0.314 

TABLE  IV 

Cases  Considered  (see  See.  10.2) 


Incompressible  Flow 

Compressible  Flow 

Rigid  aircraft 

Flexible  wing  only 

Flexible  wing  and 
fuselage 

Rigid  aircraft 

Flexible  wing  only 

Flexible  wing  and 
fuselage 
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SUPERPOSITION  METHOD 


RIBID  WIN6  x 
ELASTIC  WINS  o 


EXACT  METHOD 


RIBID  WINS  + 
ELASTIC  WINS  a 


Cf  c.  2  ty  FLI8HT  ALTITUDE  10,000m  (q  • 1210  k0/m2  ) 

M  •  0.8 


Fig.  1  Load  distributions  for  eoaprssslbl#  and  incompressible  flow  detemlnsd  by  tbs 
superposition  Mtbod  and  by  tbs  exact  aetbod 


2' 


■ SUPERPOSITION  METHOD 


RIGID  WING  x 
ELASTIC  WING  o 


- EXACT  METHOD 


RIGID  WING  + 
ELASTIC  WING  o 


FLIGHT  ALTITUDE  10,000m  (q-1210  k0/m2) 


M*  O.G 


Pig.  2  Pitching  no— nt  distributions  for  eoaprculblc  and  incoaprsssibls  flo*  deter 
■ined  bar  the  superposition  eethod  and  by  the  exact  eethod 
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SUPERPOSITION  METHOD 


Pi*. 3  Superposition  of  angle-of*»tt»ck  distributions  (Og)  and  comparison  of 
results  with  those  of  exset  method 


Idigr—s) 


&t 

(dfffrtesj 


Pig.  4 


Elevator  deflections  as  a  function  of  CL 
(CMf  =  0) 
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6f 

(digrves) 


st 

(df/r—s) 


Pig.  5  Elevator  deflection  u  a  function  of  CL 

(C-  =  -0.015//(1  -  E2) 

"f 


Pig. 6  Shift  of  centre  of  preeeure  due  to  eeroelastlcity,  with  and  without 

coapreeeibilltjr  effect* 


0  02  0.4  0.6  0.9  M  1.0 

Fig.  8  Influence  of  compressibility  end  aeroelssticity  on  wing  lift  gradient 
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Sg 


(degrees) 
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hg 


(degrees) 


Pig. 10  Elevator  travel  per  ‘g* 
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